Superconducting proximity effect in interacting double-dot systems 
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We study subgap transport from a superconductor through a double quantum dot with large on- 
site Coulomb repulsion to two normal leads. Non-local superconducting correlations in the double 
dot are induced by the proximity to the superconducting lead, detectable in non-local Andreev 
transport that splits Cooper pairs in locally separated, spin-entangled electrons. We find that the 
I—V characteristics are strongly asymmetric: for a large bias voltage of certain polarity, transport 
is blocked by populating the double dot with states whose spin symmetry is incompatible with the 
superconductor. Furthermore, by tuning gate voltages one has access to splitting of the Andreev 
excitation energies, which is visible in the differential conductance. 

PACS numbers: 73.23.Hk, 74.45.+c,03.67.Bg 
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I. INTRODUCTION 

Advances in nanofabrication techniques allow nowa- 
days the realization of hybrid nanoscale structures made 
up of quantum dots tunnel coupled to normal and super- 
conducting leads r^— which provides the possibility to 
study superconducting correlations in quantum-confined 
systems with few interacting degrees of freedom in the 
presence of non-equilibrium. Subgap transport in these 
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systems is sustained by Andreev reflection 
cess consisting of the transfer of two electrons with op- 
posite spin from a normal system to the superconducting 
condensate, or conversely, the transfer of the two elec- 
trons in a Cooper pair from the superconductor to the 
normal region. Since Cooper pairs in a BCS supercon- 
ductor are spin singlets, the latter possibility has been 
proposed as a way to inject entangled electrons in a nor- 
mal system^ Furthermore, if the two electrons stem- 
ming from a Cooper pair are separated spatially, for 
example because they end up in different leads, non- 
local correlations, similar to those discussed by Einstein, 
Podolsky and Rosen^ are generated. This type of non- 
local Andreev reflection is often referred to as Crossed 
Andreev Reflection (CAR).— ~— Andreev reflection in 
quantum dots has been the object of several theoretical 
studies.^ii 

The use of double-dot systems for studying supercon- 
ducting correlations in quantum dots introduces several 
advantages over single-dot devices. While the (equi- 
librium) proximity effect is suppressed by the onsite 
charging energy in single-dot devices, non-local correla- 
tions between spatially separated dots are more easily 
achieved. For this reason, double-dot systems are ideal 
to study CAR, as the strong on-site Coulomb repulsion 
in the dots suppress the probability of two electrons en- 
tering the same dot and therefore local Andreev reflec- 
tion. As a consequence, an appropriately biased double 
dot acts as a Cooper-pair splitter. A second advantage 
is the enhanced control on the Andreev excitation spec- 



trum. In particular, a splitting of the Andreev excitation 
energies can be obtained by means of gate voltages only 
and detected by a differential-conductance measurement. 
Recently, two different versions of the double-dot Cooper 
pair splitter have been realized experimentally ) 1 1 1 1 2 show- 
ing evidence for crossed Andreev reflections tunable via 
gate voltages. 

In this paper, we study subgap transport properties of 
a double quantum dot tunnel-coupled to one supercon- 
ducting (S) and to two normal leads (L, R) as shown in 
Fig. [TJ We focus on the regime of large superconduct- 
ing gap and dots' on-site Coulomb repulsion, where only 
non-local Andreev reflection is allowed. The transport 
characteristics strongly depend on the polarity of the ap- 
plied bias voltage. The transfer of Cooper pairs from the 
superconductor as entangled electrons into the normal 
leads occurs via a singlet state in the double dot. When 
reversing the bias voltage, however, triplet states may be 
populated in the double dot. Those have a spin symme- 
try that is incompatible with the superconducting lead, 
and transport is blocked. The onset of this triplet block- 
ade phenomenon is signaled by a pronounced negative 
differential conductance. 

Superconducting correlations in the double dot can be 
described as a coherent superposition of the empty and 
doubly-occupied double-dot states. The gate voltages ap- 
plied to the individual dots allows for an independent 
control of both the detuning between empty and double 
occupation and the energy difference between the two 
states of single occupation. The latter gives rise to a 
splitting in the excitation spectrum, which in turn is de- 
tectable in the differential conductance. 

The present article is organized as follows. In Section 
ITT1 we introduce the model of the system and the formal- 
ism used to compute the Andreev current. In particular, 
we present an effective Hamiltonian for the double dot 
coupled to a superconducting lead with infinite pair po- 
tential. In Section Hn] we discuss the main results of the 
paper: the triplet blockade and the splitting of the An- 
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FIG. 1: (Color online) Schematic setup of a double quantum 
dot, tunnel coupled to two normal and one superconducting 
leads. The double-dot system can be driven out of equilibrium 
by applying a finite bias voltages to the normal leads. 

dreev excitation energy spectrum. Finally, in Section IIVI 
we give some concluding remarks. 

II. MODEL AND FORMALISM 



The superconducting lead is coupled to both dots and 
the corresponding tunneling Hamiltonian reads 

Htunn.S = ^ [ VsLC Sk<j d L,a + VsrC^cZr^ + H.C. . 
kcr 

(4) 

Finally, the Hamiltonian of the three-terminal system 
can be written as H = H ddot + H tuniltN + H tunn ,s + 

^2 n =L,R,S Hv 

We assume that the normal-state density of states in 
the leads p n is constant in the energy window relevant 
for transport. We define the following tunnel-coupling 
strengths: T NL = 2irp L \V NL \ 2 and T NR = 2irp R \ V NR | 2 
for the left and right normal leads, respectively, and 
T S L = 2irp s \V SL \ 2 and T S r = 2t: ps\V S r\ 2 for the cou- 
pling between the superconductor and the left and right 
dot, respectively. 



The system under investigation is schematically shown 
in Fig. Q] It consists of a superconductor tunnel coupled 
to two quantum dots; each dot is contacted by a tunnel 
barrier with a normal metallic electrode. We describe the 
double dot by two spin-degenerate levels with intcrdot 
and intradot Coulomb-repulsion energies: 



a=L,R 

(7,<J r 



+ 



(1) 



where a = L,R labels the left and right dot, respectively. 
The annihilation (creation) operator for an electron in 
dot a with spin a is da,a an d the corresponding num- 
ber operator n a)Cr - The intradot and interdot Coulomb 
repulsion energies are denoted, respectively, by U a and 
U. 

The- leads are labelled by the index r\ = L . R, S, which 
indicates the left and right normal leads, and the super- 
conducting lead, respectively. The Hamiltonian of lead rj 
reads 



H.c. 



ka 



(2) 

with the lead-electron operators c n ka- and c^ fc(T . Super- 
conductivity in lead S is described by means of the mean- 
field BCS Hamiltonian with a pair potential A. Without 
loss of generality, A has been chosen real and positive. 
We also take the electrochemical potential of the super- 
conductor as reference for energies and set it to zero, 
ps = 0. The tunnel coupling between the double dot and 
the normal leads is described by the tunneling Hamilto- 
nian: 



H 



tunn,7V 



ka 



[V NL > 
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Lka L 



V N RC Rk(J d R ^ + H.C. 



(3) 



A. Large-A limit 

Since we are interested in transport processes involving 
Andrccv reflection, we focus on the limit A — > oo, for 
which the superconducting lead's degrees of freedom can 
be integrated out exactly. This results in the effective 
double-dot Hamiltonian 



fffifi =-ffddot - ( d l,t d l,i + H 



a=L,R 



2 ( d R,t d l,i - d R,A + H ' c ') ' ( 5 ) 

which adds two different contributions to the double-dot 
Hamiltonian. The first one^ describes the local prox- 
imity effect on each dot and arises from local Andreev 
reflection. For the present context the more interest- 
ing contribution is the second term, which describes the 
formation of non-local superconducting correlations be- 
tween the two dots. The latter is the one which encom- 
passes Cooper pair splitting in the double-dot system. 

In realistic double-dot systems, the intradot charging 
energy, U a , amounts to several meV and is much larger 
than the other energy scales relevant for transport^ 
Therefore, we focus on the limit U a — > oo, in which 
double occupation of each individual dot is forbidden. 
This reduces the Hilbert space of the double dot to only 
nine states: the empty state |0), the four singly occu- 
pied states lever) = d^lO), and four doubly occupied 
states with one electron in each dot. A convenient ba- 
sis for the doubly-occupied sector of the Hilbert space is: 



the singlet \S) = [d R ^-d' L ^ — d R ^d J L |0), and the 
triplet states |T0) = ^= (<4.t d Li + d R,i d L,t) 1°) and 

In the limit of a large superconducting gap, the effec- 
tive Hamiltonian Eq. ([5]) couples only the empty state 
and the singlet, but not the triplet states. This coupling 
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gives rise to the cigenstatcs 

with energies E± = 5/2 ± ca, where 5 = e L + + U 
is the detuning between the empty state and the sin- 
glet. The splitting between the |+) and the |— ) states 
is 2ea = \/S 2 + 2r|, where the effective coupling Ts is 
given by Ts = V^slTsr- The Andreev excitation ener- 
gies are defined as the excitation energies of the double- 
dot system in the absence of coupling to the normal leads: 

E A >7">7'>7 — 7"A £ /2 + j'U/2 + (7) 

with 7", 7', 7 = ±1, and Ae = el — £r- 

B. Master equation and current 



where e < is the electron charge and the current rates 
Wy take into account the number of electrons trans- 
ferred to lead r\. 

III. RESULTS 

In this section, we present results for transport in the 
double-dot system in the limit of large gap and large on- 
site charging energy. We will consider the case of equal 
coupling to the normal leads, T^l — ^nr = ^n, an d of 
symmetric bias voltages, (J.l = fJ-R = M- In particular we 
will discuss the current injected in the superconductor 
Is = -h - Ir- 

A. Triplet blockade 



We integrate out the normal lead's degrees of free- 
dom in order to obtain an effective description of the 
dot in terms of its reduced density matrix. In the limit 
^nl.nr. <C ea and taking into account only processes 
in first order in Tnl^nr, the reduced density matrix is 
diagonal in the eigenbasis of the effective Hamiltonian 
in Eq. ([5]) and we only need to consider the occupation 
probabilities P x of the eigenstates \x) of H s. In lowest- 
order in r nl.nr? the occupation probabilities obey the 
master equation 



E 

x¥=x' 



{W xx -P x > 



W x , x P x ) = 0, 



(8) 



where W x » x > is the Fermi's golden rule transition rate 
from state \' to state \" ■ 

The non- vanishing rates read 



r 



77(7, ± 



2 

^NrJ 



1± 

2e A 



E±) 
E±) 



Wria^Ta =2W ?)(Tj T0 = ^Nrjfjj- (5 ~ £77) 

W T a n a =2W T0 , V cr = r NTJ f^(S ~ e n ) 

where 77 = L,R here labels both the leads and the dots, 
and 77 denotes the normal lead other than 77, f^{oj) is the 
Fermi function of the normal lead rj with electrochemical 
potential and f~{uj) = 1 — / + (cj). 

The current in lead 77 = L, R can be computed from 
the occupation probabilities by 



= I E ns* 



(9) 



hI s /{eT N ) 
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FIG. 2: (Color online) Density plot of the current injected 
into the superconducting lead 7s in units of (eTN)/fr as a 
function of the detuning 8 = el + £r + U and of the chemical 
potential fi normal lead for degenerate dots' levels Ae = 0. 
The other parameters are: Ts = 0.5(7 and ksT — 0.0117. 

If the sign of the applied bias voltage is such that 
Cooper pairs leave the superconducting electrode, then 
the pairs are split through the double dot and entan- 
gled electrons are emitted into the leads. For the oppo- 
site bias voltage, two electrons in the singlet double-dot 
state can enter the superconductor. However, if the dot 
is in one of the triplet states, transport is blocked and 
no sub-gap current can flow. Therefore, we expect the 
current- voltage characteristics to be asymmetric with re- 
spect to the applied voltage, with a blockaded region for 
voltages favoring double occupation of the double dot. 
This type of blockade is purely due to the triplet states 
not matching the symmetry of the pair amplitude in the 
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FIG. 3: (Color online) Density plot of the differential conduc- 
tance in units of Gt un n = (e 2 / h)T n / (^ivksT) as a function of 
the detuning 5 = El + Sr + U and of the chemical potential 
fi of the normal leads for degenerate dots' levels Ae = 0. The 
other parameters are: Fs = 0.5(7 and fcsT = 0.01J7. 



3. 



tunn 




-15 



S/U 



FIG. 4: (Color online) Density plot of the differential conduc- 
tance in units of Gtunn = (e 2 / h)T n / finksT) as a function of 
the detuning 8 — £l+£r+U and of the chemical potential fj, of 
the normal leads for non-degenerate dots' levels Ae = 0.25(7. 
The other parameters are: Fs = 0.5(7 and ksT — 0.01(7. 



BCS superconductor, and we refer to it as triplet block- 
ade. In Fig. f2J we show the current in S as a function 
of the normal leads' chemical potential \i and of the de- 
tuning 5 for the case of degenerate level positions in the 
two dots, sl = sr- The blockaded region is apparent 
and it is denned by \i > (e L + e R )/2 + U = (S + U)/2. 
At the onset of the triplet blockade, the differential con- 
ductance, Gdiff = edls/dfi, becomes negative, see Fig- EI 
The triplet blockade is not lifted by removing the degen- 
eracy of the dot levels, since this will not introduce any 
coupling between the triplet states and the singlet. In 
the n < region of Fig. [5J Cooper-pair splitting occurs 
with a strong a peak around zero detuning due to strong 
proximity effect on the dot. The peaks of positive dif- 
ferential conductance in Fig. [3] map the dispersion as a 
function of detuning of the Andreev excitation energies. 
Since the limit of infinite on-site repulsion has been con- 
sidered, it is obvious that transport is only due to crossed 
Andreev reflection and that pairs are split with unitary 
efficiency. 



B. Andreev excitation energy splitting 

In the case that the dot level positions are non degen- 
erate, Ae ^ 0, the Andreev excitation energies are split 
as described by Eq. 0. This splitting is controlled only 
by gate voltages, but it is analogous to the Zeeman split- 
ting of the Andreev excitation energies in a single dot. 
The splitting of the excitation energies appears in the 
resonances in the differential conductance of the system, 



as can be seen in Fig. [4] where the differential conduc- 
tance is plotted as a function of the chemical potential /j, 
and the detuning 5 for a fixed value of Ae ^ 0. The large 
negative differential conductance pinpoints the boundary 
of the blockaded region. It has to be noticed, that even if 
the left-right symmetry is broken due to the finite value 
of Ae, still the two currents flowing in the normal leads 
are the same, i.e. Il = Ir- 



IV. CONCLUSIONS 

We have studied subgap transport from a supercon- 
ductor through a double quantum dot with large on-site 
Coulomb repulsion to two normal leads. Non-local su- 
perconducting correlations in the double dot give rise to 
non-local Andreev transport splitting Cooper pairs in lo- 
cally separated, spin-entangled electrons. When the po- 
larity and magnitude of the applied bias voltage is such 
that electrons entering from the normal leads populate 
a triplet state in the double dot, then the double dot 
is effectively decoupled from the superconductor, and 
transport is blocked. This triplet blockade shows up as 
an asymmetry of the I—V characteristics. By indepen- 
dently changing the level positions of the two dots, a 
splitting of the Andreev excitation energy spectrum can 
be achieved. This splitting can be experimentally ob- 
served in differential-conductance measurements. 
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